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Abstract 

The distribution of the characteristic polynomial Z{U, 6) oi N x N matrices U in the Cir- 
cular Unitary Ensemble is studied by the method of second quantization for one-dimensional 
fermions. For infinite A'^ the Gaussian distribution of Z{U, 6) is established straightforwardly 
by bosonization. A general expression for the n-point correlation function of the characteristic 
polynomial at different points is given by this method. The case of finite A'^ is discussed. 

The statistical properties of the Riemann zeta function Q have been extensively studied analyti- 
cally 1^ and numerically 0| and their analogy with the corresponding properties of ensembles of 
random matrices was investigated within the framework of the random matrix theory |^ . Recently, 
the distribution of values taken by the characteristic polynomials 

N 



Z{U,e) ^dei [I -Ue'") = W h -e*(''-^^n (1) 



of N X N unitary random matrices U with eigenvalues e~*^^ belonging to the circular unitary 
ensemble (CUE) was investigated in In particular, it was shown by explicit calculations 

that in the limit N oo the distribution of real and imaginary part of log Z{U, 6) divided by a 
factor ^(1/2) log AT is a standard normal distribution in two dimensions. This was conjectured to 
mimic the similar behavior of the Riemann zeta function high up the critical line. The convergence 
of the corresponding cumulants to the gaussian limit as A — s- oo was also explicitly calculated and 
they were conjectured to describe the corresponding properties of the Riemann zeta function. 

In this paper I use the equivalence between the CUE and a theory of fermions in one dimension 
to calculate the statistical properties of log Z{U, 9) using the method of second quantization. One 
of the most prominent features of this approach is that the case N — oo can be studied from 
the very beginning, thus avoiding the tedious finite A calculations and the asymptotic expansion. 
For infinite A the calculations simplify a lot, since in this case the fcrmionic theory is equivalent 
to a theory of free bosons, the fact known under the name of bosonization In what follows I 
calculate the distribution functional of the characteristic polynomial (|l|) for infinite A using this 
equivalence. In mathematical literature the bosonization is known under the name of Frobenius 
formula for irreducible characters of the permutation group or Szego asymptotic formula for 
Toeplitz determinants. 

Let me briefly describe the relation between the CUE and the fermions in one dimension. In 
the random matrix theory one is interested mainly in calculating statistical averages of functions, 
which depend only on the eigenvalues e"'^^ of U. Consider some symmetric function of eigenvalues 
f{9) = f{Oi, . . . , 9n), and its average, defined as 

(/)-(^^/^'^^ \D{e^')\'f{e), (2) 
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where the (Haar) measure of integration is defined with help of the Vandermonde determinant 

Diz) ^ Diz,, . . . , zat) = det [#"1^. ,^, 2 = Y[{z, - z,) (3) 

j<k 

Consider a quantum particle on a ring < 6 < 27r, described by the wave-function of the n-th 
orbital: 

MO) = 4=e"'- (4) 
y Zn 

The Vandermonde determinant D (e*^) is therefore proportional to the Slater determinant: 

m ^ MOu . . . , = det (e'(--)«.) (5) 

composed of particles (fermions) occupying the orbitals n = N — k = 0,...,N — 1. The propor- 
tionality factor coincides exactly with the square root of the normalization factor in front of the 
integral in (||) and this average can be rewritten as a quantum-mechanical expectation value: 

(/> - (*ol/ l*o> = J d'^e me) fie) mo) (6) 

of the operator / defined as f{6) in the coordinate representation. This correspondence of the RMT 
and one-dimensional fermions will be used extensively throughout the paper and in particular the 
statistical average (. . .) and the quantum expectation value (^o| • • • I'^o) will be interchanged in 
the course of the paper by the virtue of (|^). Some application of the fermionic picture will be 
presented in what follows. 

The central object of my discussion is the logarithm of the characteristic polynomial (|l|): 

N 



m = log ZiU, 0) = ^ log 1 - e^e^-*^^) . (7) 



Expanding the logarithm, the equation (0) can be rewritten as 



fc=l j=l k=l 

where we have used the Fourier transform of the density operator: 

TV „27r N +00 

p, = tr {U") = ^ e-^''^ ^J'^^ e-^''pie), p{d) ^T.^iO-e,) = — ^ e'^^Vfe. (9) 

j — l ^ j=l k— — oo 

In order to calculate the statistical properties of L{9) using the correspondence between sta- 
tistical average and quantum expectation value it is convenient to employ the method of second 
quantization. We introduce creation and annihilation operators cj^ and c„ for a fermion on the 
n-th orbital with usual anti-commutation relations: 

The quantum state j^I'o) is then defined by the action of the creation operators on the vacumir: 

N-l 

l^o) - n ^«|Vac>. (11) 

n=0 

The second-quantized form of the density operator pk is given by the standard rules [lOt] : 

+00 



n— — oo 
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for fc ^ 0, while for fc = the definition (g) gives po = N — the total number of particles. When 
acting on the state |^'o) the density operator pk creates a linear combination of states in which 
one particle on the n-th orbital is moved k orbitals down, provided the orbital n — k is empty. 
This simple observation yields the important result (for k ^ 0): 

{pkpD = 4,pC2(fc) - 6,,,{^o\Pkpi\^o) = 4,p X I 1^^' j^j J ^ (13) 

otherwise obtained by using the Wick theorem. The correlation function C2(fc) = {pkpl) with 
(^2(0) = N"^, is nothing but a Fourier transform of DOS-DOS correlation function Q of CUE : 

1 +°° AT / /V \ 2 / sin M \ 2 

C,{e;N)^{p{e)pHO))^- J2 e''''C2ik;N)^-6{e)+(-) - i —-^ ] . (14) 

The symmetry of the correlation function C2(fc; A^) with respect to fc — s- —k follows from the fact 
that for finite N the density operators pk and commute. For infinite N , as we shall see, this is 
not true. 

I now proceed to calculate the correlation functions of the logarithm of the characteristic 
polynomial. In the work |^ the correlation functions of L(Q) were calculated at the same point. 
Here I generalize this result, to the case of the two-point correlation function, and calculate 

A'2(0i-e2;iV) = (£(0i)Lt(02)) (15) 
P^(er-d2\N) - (L(0i)i(02)). (16) 

where in the coordinate representation L'^iff) — L*{9). 

Due to the translational symmetry all the correlation functions depend on the difference 9 — 
9i — 92 only. From (15) and ( p^ the correlation functions of real and imaginary part of L{9) can 
be easily obtained. I notice that P2{9) vanishes due to the fact that the Kronecker delta in ( p^ ) 
is never satisfied. Moreover, it follows that 

Q2{9^~92;N) ^ (Re£(gi)ReL(g2)) ^ ^""^^^^^^ ^^'^^ 

i?2(0i-02;iV) ^ (Im£(gi)Im£(g2)) = ^^-^^(^i^- ^2; A^) ^^^^ 

and in addition there exists a cross-function given by 

X2(gi-g2;iV)= (Im£(0i)ReL(g2)) =-(Re£(gi)Im£(g2)) = ^"'^^^^^^^ ^^'^^ (18) 

It can be checked that {L{9)) = (L^{9)) = so there is no difference between connected and 
disconnected correlation functions. 

The calculation of K2{9;N) using the definition (^) is straightforward: 

°° ike °° ike N ike °° ike 

K2i9,-92;N)= Y: ^W]i> = E^^2(^;^)-E — + ^ E ^ (19) 

k,p=l ^ k=l k=l k=N+l 

For correlations at the same point, 9 — 0, the result is real and its half coincides precisely with 
the expression (43) in 

k=l k=N+l 

which behaves as (1/2) log A'' for large N. It is worth noticing that for 9^0 the function K2{9) is 
complex, therefore there exists correlation between real and imaginary parts of L{9) at different 
points, a fact which is missed when the correlation function is calculated at the same point. 
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Now the calculations for the whole distribution of L{9) and L^{d), equivalent to the distribution 
of TieL{9) and Im L{9) is presented. In order to be able to calculate general n-point correlation 
functions, one has to consider the following generating functional: 



(exp ^ [s*i9)m + si9)L\9)] |^ (21) 



where s(9) and s*{9) are the source terms. Any n-point correlation function can be represented 
as a functional derivative of S[s, s*]: 



{m)m) . . . L*i9„,)L*i9„.+,) . . . E[s, s* 



(22) 

s,s*=0 



where Sm, stand for s{9m) and s*{9i) respectively. Using the Fourier transform of these source 
terms 



2tt 

ike , 



Sk= / d9e-''''s{9) , (23) 
Jo 

and the expression (|^) the generating functional (|2l| ) can be rewritten as 



E[s, s*] = (exp I - E ^ {'kPk + skpi) I j ■ (24) 

It will be calculated in the limit N = 00. It is convenient to redefine the numbering of one-particle 
orbitals so the upper occupied level in {"^q) corresponds now to n = and all the states with 
n < are occupied. This state is the infinitely deep Fermi sea — the ground state of the fermionic 
system, if one-particle energy En is an increasing function of the level index n. In addition, this 
state is now annihilated by the action of pk for fc > 0: it is impossible to promote a fermion from 
the state n to an empty lower state n — k. It is well known from the theory of one-dimensional 
correlated electrons that the operators pk and pi, acquire non-zero commutation relations in 
the presence of infinite filled Fermi sea (Schwinger terms). To show it let us begin with 

+00 +00 

[Pk,pl]^ J2 ['i-k('n,cl+pCm]= J2 (ci-fcC«-p-4_fc+j,C„) . (25) 
n,m— — oo n— — oo 

The result is nonzero, since the operators are not well-behaved. It is necessary to extract the 
singular part by introducing the normal ordered operators in the state |5'o) by extracting the 
expectation value in this state: 

• ^In^n '■— (^In'^n ~ {^In^ri) — cl^Cn — Sm,nnjn 7 (26) 

where nj^ = {clnCm) is the Fermi-Dirac distribution of occupation numbers in the ground state. 
Rewriting cl^Cn ='■ cj„c„ : +5m,nn'^ and substituting it into (|2^) the commutation relation of 
density operators becomes 



-t-00 +00 



n— — 00 



In this expression the normal ordered operators were canceled, since they are not singular. The 
last sum equals to the number of orbitals from 1 to k. 

When calculating the matrix elements as in ( p^ ) the order of pk and p^ is important when 
these operators do not commute. One requires that the matrix element should coincide with the 
statistical average in the coordinate representation (^. Suppose that an operator pk for fc > 
happens to be next to the left of the state Iv^o)- The result would be zero, since this state is 
annihilated by pk for fc > 0, which is not true for the statistical average. In order to obtain the 



4 



correct expression all operators pk for fc > must be placed to the left of p\. for fc > 0, in the 
so-called anti-normal order. The operators in the generating functional ( p4[ ) must be therefore 
rearranged as 

E|.,..l = (exp(-|:£^)exp(-|:ifi)) (28, 

Using the fact that the commutator of pk and p\ is a c-number and pk annihilates the ground state 
for fc > 0, the generating functional is given by the famous Baker-Campbell-Hausdorff formula: 

= (cxp(-^:^)expf-f:f|i))=exp(f:4|^[p.,pt]) 




(29) 



which is obviously gaussian. In particular, taking the corresponding derivatives of S[s, s*] with 
respect to Sk/k and s^/fc in analogy to ( p^ the correlation function of powers of of pk or traces 
of can be calculated: 

\fc p I \ k p Ik 

in accordance with the results of jl^. This result and the gaussian characteristic functional (|2^ ) 
are consequences of the anomalous commutation relations ( p7| ) that are fulfilled by the operator 
Pk and its hermitean conjugate and the fact that pk annihilates the ground state for N = oo. 
It is worth mentioning here that the correlation function (|o|) and its generating functional ( ]29| ) 
are in fact the restatements in terms of the field theory of the so called functional central limit 
theorem for log Z discussed in . There it was proven using the Frobenius formula for the 
irreducible characters of the permutation group, which constitutes the mathematical basis of the 
bosonization 

Returning to the function s{0) of the angle, the bilinear form in the exponent of ( p9[ ) can be 
rewritten as a double integral 

l0gS[s,S*]=/ / ^^S*(^i)i^2(^l-^2;OO)S(02), (31) 

where the "propagator" is given by (19|) for infinite iV, i.e. 



K2{e- 00) = V ^e-^'^" = - log (1 - (32) 

k=l 

and the infinitesimal imaginary part ry of the angle 9 was included to ensure convergence a,t 6 ^ Q. 
It was not necessary for finite N ^ since in ([l9| ) the first sum is finite and the second one is 
absolutely convergent. It is well known ||, |l^ that an ultraviolet cut-off l/rj should be introduced 
when dealing with an infinite Fermi sea. Using this cut-off, the correlation function ( ^2[) at the 
same point {6 = 0) is described by 

if2(0;«3) = -log(l-e-2'') «logi-, (33) 

which should be compared with the leading behavior of K2{0]N) ^ logA^. In fact, for large but 
finite A'^ one can make the scaling ansatz for the correlation function K2{0; N,r]): 

i^2(0;iV,77) = logi- X / (^log TV/log (34) 
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Figure 1: The correlation functions Q2{9;oo), i?2(6';oo) and X2{0;oo) for rj = 0.001. 



with f{x) ^ 1 as a; — > cxD and f{x) « x for x — + 0. The scahng form can be justified by 
calculating the leading behavior of (|l9|) for iV — > oo and = 2iri 0. Using the scaling ( |3^ ) 
the correspondence between TV and l/2r] is established for the leading behavior (in logA^) of the 
second moments of TieL{9) and lmL{9) at the same point. Whether the bosonization, which 
applies for the case N = oo only, can provide results for other correlation functions for finite N 
through some scaling relations like (^) is an open question. 

Finally I give explicit expressions for correlation functions Q2, R2 and X2 for 6* 7^ calculated 
using (0,|8h. Separating the real and imaginary parts in (p3) they read: 



Q2(^;oo) =i?2(^;oo) = ilogf ) (35) 

2 \^|2sm2|y 

X2{e; 00) = mod I (36) 

where the mod7r/2 means that the values taken by ^2(6'; 00) lie in the interval [— 7r/4,7r/4] so 
it is a periodic function of the angle. The functions ( p5|j3^ ) are shown in Figure |[ It would be 
interesting to compare these correlation functions with the corresponding functions for real and 
imaginary parts of the Riemann zeta function high up the critical line. The function i?2(^; 00) is 
related to the variance of number of eigenvalues in the arc of length 9 which was studied in [ pT| . 

In this paper I have shown how the one to one correspondence between the probability measure 
and Slater determinants of fermions in one dimension can be used in order to calculate different 
statistical properties of Circular Unitary Ensemble of random matrices. In particular, the cor- 
relation functions of the density of states, which corresponds to the density of the fermions are 
obtained with help of the Wick theorem as a diagrammatic expansion. This method was applied 
to the calculation of correlation functions of characteristic polynomials. The case of = 00 was 
treated by the method of bosonization and the generating functional of correlation functions was 
calculated exactly yielding an alternative derivation of some of the results of 0, |ll|, ^ . The 
corresponding distribution was shown to be gaussian corresponding to the functional central limit 
theorem of |0] . In the future it would be interesting to investigate by the present method the case 
of finite in order to obtain the correlation functions of the characteristic polynomials at the 
same point. 

In conclusion I would like to remark that the finite N calculations presented in this paper and 
in H, 0, |ll], ^ are related to the results obtained for the (static) correlation functions of strongly 
correlated particles in one dimension. Indeed, it is known that the different circular ensembles 
of random matrices correspond to the Calogero-Sutherland model , which describes a system 
of N interacting particles. This model with interactions falling off as the inverse square of the 
distance between the particles and proportional to the coupling constant A, was shown in |13] 
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to be equivalent to the system of effective free particles with fractional quantum statistics. The 
values of the coupling constant A = 1/2, 1 and 2 correspond to the orthogonal, the unitary and the 
symplectic ensemble respectively. In the present work only the A = 1 Calogero-Sutherland model 
equivalent to the free fermions was considered within the framework of the second quantization. 
The generalization of the present approach to other ensembles, corresponding to more exotic 
fractional statistics would be an interesting issue. 
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